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Abstract 

An analytical approach to quantum mechanical wave packet dynamics of laser-driven particles is 
presented. The time-dependent Schrodinger equation is solved for an electron exposed to a linearly 
polarized plane wave of arbitrary shape. The calculation goes beyond the dipole approximation, 
such that magnetic field effects like wave packet shearing are included. Analytical expressions for 
the time-dependent widths of the wave packet and its orientation are established. These allow for 
a simple understanding of the wave packet dynamics. 
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I. INTRODUCTION 



Experimental progress in laser physics has enabled to reach intensities that can accelerate 
electrons to velocities comparable to the speed of light. The theoretical description, i.e. the 
solutions to the Dirac equation describing particles that are exposed to plane light waves 
( Volkov states) have long been known 1] . Since a particle that occupies such a Volkov state 
is not localized in space, the dynamics is neither very intuitive nor does it describe a realistic 



situation. Different methods have been developed to describe 



ocalized wave packets. One 



approach is to numerically solve the Dirac equation as in P, \q |4|, |5( . A further possibility 
is to superimpose Volkov solutions to find the dynamics of wave packets f|. It has been 
found that the wave packets deform and spread non-isotropically for high laser intensities. 
Although this behavior is seen when the wave packets are plotted, there is no detailed 
understanding of their dynamics. 

Going to lower intensities where the particle dynamics can be described by means of the 
Schrodinger equation the dipole approximation can often be employed. The consequence is 
that the magnetic field is neglected. This case can easily be handled analytically (see for 
example 7f), but without the magnetic field the deformation effects do not appear and the 
wave packet spreading reduces to that of a free particle. 

The purpose of this paper is to establish the wave packet dynamics in the regime where 
the dipole approximation is not valid anymore and the velocity of the particle is still slow 
compared to the speed of light, though non- negligible. That means, the particle can be 
described by the Schrodinger equation which includes all relativistic effects of order 1/c. 
Since the magnetic field plays an important role in the regime considered here, deformation 
effects are to be expected. The wave packet will be characterized by the trajectory of 
its maximum, the principal axes and its widths. The establishment of the corresponding 
analytical expressions describing the wave packet will give intuitive insight into its dynamics. 

Having defined the laser fields and the appropriate field expansion in the first section, the 
dynamics will be solved for a classical particle. Later on, this will be used for a comparison 
with the trajectory of the quantum mechanical wave packet. In section HVI to lVll the quantum 
mechanical wave packet dynamics will be deduced. After solving the Schrodinger equation 
for a spacially restricted area, a gaussian wave packet will be constructed. In the last section 
the results will be applied to a particular example of a short laser pulse. 
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II. LASER FIELD EXPANSION 



The following calculation will refer to the dynamics of an electron, that is exposed to 
a laser pulse. The emitted light is supposed to be a linearly polarized plane wave with 
arbitrary shape. If the pulse travels in z-direction, the vector potential is given by 

A = x ■ A{ujt - kz) = x ■ A(ip), (1) 

where the polarization is in x-direction. The wave vector k and the angular frequency u are 
as usual connected by the speed of light c — uj/k. With the assumption that the particle be 
close to the plane z = 0, i.e. (kz) 2 <C 1 the vector potential can be expanded: 

A « x [A(ut) - kz ■ A' (tot)] . (2) 

The prime denotes the derivative with respect to the phase ip, whereas a dot will be used 
for time derivatives. The following example is a useful check for the validity of the approxi- 
mation: If the laser emits visible light of wave length A = 800nm, then the particle can be 
displaced up to a few hundred atomic units away from the plane z = without kz exceeding 
the order of 10 _1 . 

As in the expansion (0), the following calculations will be carried out to the first order 
in (kz). That means effects that are due to the z-dependence of the vector potential will 
be found which are not contained in the dipole approximation, where any z-dependence is 
eliminated. 

The electric and magnetic fields are derived from the vector potential: 

E = —-A = — — (A' — kz ■ A") ■ x (3a) 
B = VxA = -k-A'-y. (3b) 

As opposed to the dipole approximation, the magnetic field is not equal to zero here. 

III. CLASSICAL EQUATIONS OF MOTION 

Although the fully relativistic dynamics of a classical particle driven by a plane wave, is 
well-known (see for example jjj), a non-relativistic treatment beyond the dipole approxima- 
tion will be presented. The purpose is to find a result that can be used for a comparison of 
the trajectory of a quantum mechanical wave packet and the corresponding classical solution. 
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In the following, atomic units will be used, i.e. Planck's constant h, the electron mass m 
and the unit charge e will be set equal to one. 

The motion of a classical, non-relativistic charged particle is given by the equation 



'i=E + -xB. (4) 



c 

Applied to the fields (jSJ) one arrives at the following component form: 

x = ~(A'-kz-A") + *-z-A' (5a) 

y = o (5b) 

z = -— ■ x ■ A' (5c) 
<r 

Using to ■ A' = A equation (p)aj) can be integrated and inserted into (JST-j) : 

x = — A + — ■ z ■ A +p x (6a) 
c c 2 

z = ^(-A-kz-A' - Px )A' (6b) 
& c c 

Here the initial conditions are such that x(t = 0) = p and x(t = 0) = xq. Now, equation 
f!6b|) is formally integrated and one finds the following integral equation: 

- = ^-A 2 --^ + P --- A f(kz)A>>dt (7) 



c 2c 4 c 3 c c 4 . 
At this point, some simplifications can be made. A/c 2 and p x /c have to be a small 
number such that the particle dynamics is non-relativistic. This is seen from (J6aJ) (x/c = 
-A/c 2 +p x /c+ ■ ■■ < 1). 

Now, the maximal values of A(t) and A'(t) are assumed to be of the same order of 
magnitude which is justified by the following example. For a typical laser pulse of sinusoidal 
form with an envelope R(t) one finds 

A{t) = R(t) ■ sin(cjt) (8a) 
A'(t) = ■ sin(ujt) + R(t) ■ cos(ujt) . (8b) 

A laser pulse which is not too short will reach its maximum R only after several quarters 
of a period T, i.e. the slope of the envelope will typically fulfill the condition R(t) 
Rq/(T/4) ujRq. With that it follows from (jHJ) that the maximal values of A(t) and A'{t) 
are of the same order of magnitude. 



In the following, only those terms are kept that are not higher than second order in the 
small quanities A/c 2 , A'/c 2 , p x /c and kz. Thus the integral in (0) vanishes and the equation 
can be integrated to give an explicite solution for the motion in z-direction: 

' ' A 2 dt - ^| / Adt + p z ■ t + z . (9) 



2c 3 J c 2 _ 

This can be inserted into (jfjajl . and after neglecting the high order terms again, this can 
be integrated: 

x = ~{ 1 + ^f) J Adt+^( Pz -t + z )A + p x -t + x . (10) 
The dynamics in y-direction which is defined by (|5b|) is simply a free motion: 



y = p y -t + y . (11) 
The equations of motion found here could also be obtained as a limit of the fully rela- 



tivistic solution 



with the corresponding approximations. 



IV. QUANTUM MECHANICAL DESCRIPTION 



The quantum mechanical dynamics of an electron in a laser field is determined by the 
Schrodinger equation: 

2 -, 



9 I 



-iV - -A 

c 



_V 2 + —(A -kz-A')^- + \a 2 - -„AA' 



c 



dx c 2 



i>. (12) 



Note here that the vector potential is assumed to be in the Coulomb gauge, i.e. V • A = 0. 
The particle is supposed to be localized in the vicinity of z = 0, i.e. terms of order (kz) 2 
are neglected. This will also be employed to find the solution of the Schrodinger equation. 
It will be solved separately for the orders (kz)° and (kz) 1 using the following ansatz: 



V> = expi[p x (x - x ) +p y (y - yo) +p z (z - z ) + u(t) ■ kz + w(t)] . 



(13) 



With that the partial differential equation reduces to two ordinary differential equations, 
one for each order of (kz): 



u 



-P^ A > + La A 1 



w = — - 



2„,2 



Px + Pv + Pz + ^Pz ■ U + k U 



^■a + \a 2 ] 

c c z 



(14a) 
(14b) 
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These equations can be solved by mere integration: 



u = -^A + -^-A 2 (15a) 
uc 2ujc 2 

w = - 1 -[pl+p 2 y +pl}t + ^(l + ^) j Adt (15b) 
_ _L (i + El + i) / A 2 rft+^- / A*dt-— R [ A 4 dt. 

2c 2 v c c 2 y y 2c 5 y 8c 6 y 

Now, those terms will be dropped that are by a second order factor of A/c 2 less than other 
terms. Thus the solution of the Schrodinger equation reads 



ip = exp i 



p x (x - x ) + p y (y - y ) + p z (z - z Q ) + (^~^ A + ^~2^ 2 ) ' kz 



dt 



• (16) 



This wave function was found under the assumption that the particle be localized in z- 
direction, i.e. a wave-packet has to be constructed out of these non-localized wave functions 
to get a consistent solution. That will be subject of the next chapters. 

V. GAUSSIAN WAVE PACKETS 

A convenient method is to work with gaussian wave packets. They possess well-defined 
widths in both momentum and coordinate space, and they are rather easy to handle. 
Now, functions of the general form 

= -4=exp [if(x,p,t)] (17) 
y In 



are being superimposed with a gaussian shaped weighting factor 

1 {p-i 
l ■ exp — 



c 1 (P-Po) 2 nc x 
f = ; • ex P ( 18 ) 



to define the wave function \1>: 



*(x,t) = (TTv^Apf 172 J°° exp- (P A J 2 o)2 - exp [if(x,p,t)]dp. (19) 

Here f(x,p,t) is a general real function that depends on time, one momentum and one 
spacial variable. The maximum of the momentum distribution is po and its width is given 
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by Ap. If both functions are normalized according to 

<f)*(x,p,t)<f)(x,p,t)dx — 8{p — p) (20a) 

/oo 
C(p)aP)dp = l- (20b) 
-oo 

then this holds for \1/ as well: 

/oo poo 
^*(x)^(x)dx = / (j)*(x,p,t)C{p)(i>(x,P,t)£{p)dpdpdx = 1. (21) 
-oo J — oo 

Now the function f(x,p,t) is expanded with respect to p at the point po: 

/(p) = /(Po) + /'(po)(p - Po) + §/"(Po)(p - Po) 2 • (22) 

This can be a good approximation if the momentum p is close to po due to a narrow width 
Ap in momentum space. If f(x,p,t) is only quadratic in p, as will be the case here, the 
expansion will even be exact. Now the integral (fl"9|) can be done by means of 

/oo r^r jj2 

exp(— ax 2 + j3x)dx = y — exp — . (23) 

The wave function \l/(x) and the probability density |\I/(x)| 2 are found to be 

* w = H " ^' f ' M ) "' /2 exp (i/(po)) • exp " i -flpy k) (24a) 



It can be read off, that the maximum of the wave packet is determined by the derivative 
f'(x,p,t) whereas the second derivative f"(x,p,t) is responsible for the wave packet spread. 

VI. WAVE PACKET DYNAMICS 

The results of the last section can now be applied to construct a localized wave packet out 
of the wave functions (|16j). The formula for the probability density (|24j) has been derived 
for one dimension. However, the three dimensional problem at hand can be reduced to three 
one dimensional ones, i.e. the wave function (fTH|) can be rewritten in product form 

^{Px,P y ,Pz) = exp(p a .) • exp(p^) • exp(p 2 ) . (25) 
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For that, the mixed term p x p z in (fTp^l needs to be removed, which is accomplished by means 
of a coordinate transformation. The p^-dependence separates from the beginning, i.e. the 
transformation is only in x- and z-direction. The terms E q in the exponent (J 16)) that are 
quadratic in p x and p z are written in the following form 




where T denotes the transpose of a matrix. The 2 x 2-matrix in ()26|) possesses the following 
eigenvectors vi/ 2 with the eigenvalues Ai/ 2 : 



Vl = ^-(]) Ai = t-4 / Aft 



c 2 



(27) 

From that one finds the following transformation law: 

1 V2 Y (28) 

z = -j=(-x + z) z = —(x + z) 

The new coordinate system is related to the old one by a rotation of 7r/4 to the left. In the 
new coordinates the wave function ()16|) can be rewritten in product form: 



ip = exp % [p y (y - y ) - \p\t - J A 2 dt] 



x exp i 
x exp i 



p x {x-x -^A(x + z) + ^-jAdt-^J A 2 dt) -lpl(t-^J Adt) 
p z (z-3o + + z)-7Tcf Adt -^J A2dt ) ~ 5# (* + h I Adt ) 



(29) 
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Applying (j24b|) to this wave function one finds the following expression for the probability 
density: 



x 



x 



(x,y,z 



V2^ 



i + 



-1/2 



Ap; 



2 



exp 

-1/2 



L [y - yo - p y t] ' 



Ap 2 



1-1/2 



M 

2 



exp 



x 



x - ^A(x + z) + ^jAdt-^f A 2 dt -p x (t-±j Adt) (30) 



x exp 



M 

2 



(t-±fAdt)) 

z + ^A(x + z) - fef Adt- / A*dt-p z (t+j;fAdt) 



(t+ifAdt) 



The motion of the maximum x m (t) of this gaussian wave packet is easily found by equating 
the arguments of the exponential functions to zero. It is instructive to transfer these equa- 
tions back to the original coordinates by means of (|28jl. The following equations are found: 



x 7 



-in 

c 



Pz 

c 



Adt + — ■ z m ■ A + p x ■ t + x 

c 2 



S Adt + ' 1 + ^ ' A + Px ' 1 + Xo 

Adt + p z -t + z 



2c. ' A dt 

ym=Py-t + y . 



P, 



(31a) 

(31b) 
(31c) 



These are identical to the equations of motion of a classical particle ()9I10I11|) . 

To fully understand the wave packet dynamics it is helpful to examine the directions of 
the principal axes (i.e. the directions of those axes where the widths of the wave packet are 
extremal) and the widths of the gaussians along these axes. The probability density (J3*U|) 
relative to its maximum is proportional to 



\ip\ 2 oc exp — 



&[x-£,A(x + z)]' ^[z + ^Mx + z 



Z " _| Z L AC" v 'A _| Z L \ 1 J 

l + (^(t-±fAdt)) 2 l+(^(t + ±jAdt)) 



(32) 



The coordinate transformation needed to find the normal form of the exponent will only con- 
cern the x- and the z-direction. Although the exact transformation can be found the result 



is complicated and not very enlightening. It is much more instructive and less cumbersome 
to consider an approximative transformation for the case that the initial wave packet is 
axially symmetric (with respect to the y-axis), i.e. Ap x = Ap z = Ap. Furthermore, the 
terms 1/c 2 J Adt are considered to be small compared to t. This is a reasonable assumption 
for two reasons. Firstly, the dynamics is non-relativistic, i.e. (A/c 2 ) 2 <C 1 (see section ITTT]) 
and secondly, the integral partly cancels if the vector potential describes an oscillating laser 
pulse. With that, the following quadratic form QF needs to be transformed to its normal 
form: 



QF 



Ap 2 
2 



1 + 



Ap- 



i V / i x2 



(33) 



Now, the following orthogonal transformation is applied: 

To the order (A/c 2 ) 1 this is a rotation of the coordinate system by an angle of A{t)/(Ac 2 ) 
to the left. With (A/c 2 ) 2 <C 1, the quadratic form (|3*3*jl is diagonalized and the probability 
density becomes 



oc exp 



Ap?, 



•y 



Ap 2 
2 



Ap2_, 
2 L 



1 



A) X + [1 



1 + m 1 + 

The corresponding widths (for the definition see (jlHjl ) are given by 

2 



A Z 



■ (35) 



AX' 



AZ 2 



Ay 2 



Ap 2 

2 
Af? 

2 



( 1+ i- 4 ) 










Ap 2 



With that the ratio of AX and AZ is found to be 

AX 1 „ 



(36a) 
(36b) 
(36c) 

(37) 



VII. DISCUSSION 



The equations of motion for the maximum of the laser-driven wave packet ([31)1 . which are 
identical to the classical motion, contain several terms that are proportional to the vector 
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potential A(t) or its integral J A(t)dt. Typically, these are terms that oscillate with the 
laser frequency lu, such as the term responsible for the motion in E-field direction (-1/c 
J A(t)dt). A different term J A(t) 2 dt is contained in the equation for the z-direction, which 
is responsible for the particle drift in laser-propagation direction. Since the integral is over 
the square of the vector potential this term grows as the laser pulse passes, no matter what 
shape it has. That means, a particle will always be displaced in propagation direction of a 
laser pulse. This is understood by the magnetic Lorentz force that is excerted on the particle 
as it moves in polarization direction. It is necessary to go beyond the dipole approximation 
to find this term. Finally, there is a constant motion term (p ■ t + const) for each direction 
that is determined by the initial conditions. 

The wave packet itself is characterized by its principle axes. Their directions are those 
of the basis vectors of the initial coordinate system which is turned by an angle of 7r/4 + 
A(t)/{Ac 2 ) about the y-axis. The widths of the gaussians along the principle axes are given 
by (jSEJ). The terms containing the vector potential are only (typically oscillatory) corrections 
to the wave packet spread Ax of a free particle 




Ap 2 \ V 2 

which is determined by the initial momentum width Ap. Thus the wave packet shows the 
well-known behavior, that the narrower the packet is initially, the faster it spreads due to 
the corresponding greater width in momentum space. 

In the following example a gaussian shaped laser pulse of a few cycles is considered 
(FIG.HJ). It is given by 



A(t) = Rc 2 sin(u;t — ip ) exp 



Ldt — Cfo X 



# = 0.25, Az = 6, ifi = 3.5- Az (39) 



Az 

The total area of that anti-symmetric pulse is zero, i.e. Adt = 0. Thus according to 
(jDpiOpilj) the position of the particle after the pulse has passed is given by 

x = p x ■ t + xo (40a) 

z= — [ A 2 dt + p z -t + z ^ 282.4 + p z ■ t + z (40b) 
2c d J 

y = p y -t + y , (40c) 

i.e. it is travelling at its initial velocity and it is displaced by roughly 280 a.u. in z-direction. 
The electron does not pick up any momentum, because the final and initial values of the 
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FIG. 1: The vector-potential in the dimensionless form A/c 2 , which represents the laser pulse 
shape. Except for a factor of four, this curve also shows the time-dependent angle between the 
principle axes of the electron wave packet (in the x-z-plane) and the diagonals of the coordinate 
system (x = ±z). The angle is given by Aj(\c 2 ). 

vector potential are identical. For the same reason, the widths of the wave packet ()36|) 
reduce to the form of a freely travelling particle (jH8j) . 

The dynamics of the particle in x- and z-direction is depicted in (FIG. EJ). (A) shows 
the motion of the maximum of the wave packet. Now, the initial conditions are chosen 
such that the particle rests at the origin at the beginning. The widths in momentum space 
are Ap x = Ap z = 0.05 a.u.. The wave packet is shown in (B) for the times when the 
maximum is located at the positions marked in (A). These correspond to the times tj = 0, 
wtu/vi — ifo = ±7r, ojtm/ v — ipo = ±7r/2 and uti V = tpo, where <fo determines the maximum of 
the envelope (see It is seen that the wave packet, while spreading, deforms according 

to (|37j) . 

We note finally that pure quantum effects due to self interferences with negative Wigner 
functions are possible for the described scenario 0. 
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FIG. 2: Electron wave packet dynamics: (A) shows the trajectory of the maximum of the wave 
packet. The labels (i) to (vi) mark the positions for which the wave packets are plotted in figure 
(B). The sides of each of the squares have a width of 70 a.u. In the pictures (iii) and (v), the 
principle axes of the wave packets are shown. The directions deviate from the diagonals by an 
angle of A(t Ui / v ) / '(4c 2 ). The other wave packets in (B) are axially symmetric. 
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